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Abstract 



Vacuum asymptotically flat Robinson-Trautman spacetimes are a well 
known class of spacetimes exhibiting outgoing gravitational radiation. In this 
, paper we describe a method of locating the past apparent horizon in these 

■ spacetimes, and discuss the properties of the horizon. We show that the past 
' apparent horizon is non-timelike, and that its surface area is a decreasing func- 

Q\ , tion of the retarded time. A numerical simulation of the apparent horizon is 

""^ ■ also discussed. 

cr 

6JT)! 1. Introduction 

Vacuum asymptotically flat Robinson-Trautman spacetipes have been the subject 
r> | of much study since their discovery over thirty years ago.113 They possess some very 

■ nice features which make them amenable to analysis. They are in some sense the 
simplest asymptotically flat solutions which exhibit gravitational radiation, albeit of 
a fairly specialised nature. The feature of the spacetime which simplifies analysis is 
that the full spacetime can be built up from the solution of a fourth order parabolic 
equation on a 2+1 dimensional manifold. This is a consequence of the fact that the 
coordinate system uses a retarded time coordinate, so that in effect initial data is 
prescribed on a null hypersurface. 

Early studies of the Robinson-Trautman equation related to behaviour of solutions 
of the linearised equation. lo Lukacs et.altJ studied the behaviour of solutions of 
the full nonlinear equation, using the concept of Lyapunov stability to establish that 
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global solutions, if they existed, would converge to the static Schwarzschild equilib- 
rium. Ajjaumbcr of authors subsequently focused on the existence of such solutions. 
Schmidlllj showed local existence of solutions for sufficiently differentiable but oth- 
erwise arbitrary initial data. Rendallia showed global existence for sufficiently small 
initial data, antipodally symmetric on the sphere, and his proof was extended bjj 
SingletorCJ to remove the requirement of antipodal symmetry. Finally ChruscielB 
proved semi-global existence of solutions, for arbitrary smooth initial data. Numer- 
ical studies by Perjescj and SingletonRil demonstrate vividly the evolution of an 
initially perturbed spacetime to the steady state. It is interesting to note that the 
evolution equations exhibit exponential divergence from arbitrary initial data in the 
"backwards" time direction, and that no solutions exist in this direction.EI It should 
be emphasised these results, and the subsequent discussion, only apply to asymp- 
totically flat RobinsonpTrautman spacetimes with regular topological S 2 surfaces: 
Perjes and HoenselaersH have shown that there exists a class of Robinson- Trautman 
spacetimes with cusp singularities which evolve to static C-metrics. 

Thus the evolution of the Robinson- Trautman spacetimes would seem to be thor- 
oughly understood. We would expect then that the behaviour of physical features 
of the spacetimes should not be too difficult to pin down. Near J + , the spacetimes, 
because of their algebraically degenerate structure, could reasonably be interpreted 
as describing purely outgoing radiation around a black hole source: e.g. the decaying 
tail of the radiation after the black hole has formed (Note that there is no ingoing 
radiation, which means there can be no backscattering or similar self-interaction 
of the radiation field). However, in this paper, we will be studying the structure 
of past apparent horizons in vacuum Robinson- Trautman spacetimes; thus we will 
focus our attention on the behaviour of the "white hole" region of the spacetimes. 
Since these spacetimes evolve to the Schwarzschild geometry, the future apparent 
horizons coincide with the future event horizons at u = oo (r = 2m in Schwarschild 
coordinates). 

Our study of the apperent horizon structure of these spacetimes was initially mo- 
tivated by an investigation of the stability of the related class of electrovac spac£i 
times, the Robinson- Trautman Einstein-Maxwell spacetimes, discussed elsewhere Ed 
The asymptotically flat vacuum Robinson- Trautman spacetimes turn out to have 
a particularly well behaved apparent horizon which illustrates several general the- 
orems concerning apparent horizons. There are in fact very few well understood 
exact solutions which can be used to illustrate such theorems - there being few 
examples of non-static black hole type solutions. The asymptotically flat Robinson- 
Trautman spacetimes, in some sense the simplest asymptotically flat spacetimes 
admitting gravitational radiation, should be a useful example for the study of ra- 
diating black hole spacetimes in general and many of the general theorems relating 
to non-stationary spacetimes. The basic equations and features of the spacetimes 
arc outlined in Section 2. Section 3 contains the results concerning the past appar- 
ent horizon and its properties. Some numerical demonstration of these results is 
presented in Section 4. 
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2. Robinson- Trautman vacuum spacetimes 



2.1. Basic equations and notation 



The line element for the Robinson- Trautman class of spacetimes is given by: 

ds 2 = 2Hdu 2 + 2dudr ~ (1) 

where H = — r(lnP) m — ^, and K = A(lnP). The vacuum Robinson- 

Trautman spacetimes are foliated by a two parameter family of 2-surfaces S u ^ r 
which, for the asymptotically flat case, have spherical topology. The operator A = 
2P 2 d^ is the Laplacian on these two-surfaces. For the vacuum spacetimes, m — 
m{u) and we can use a coordinate transformation on u to make m constant. We are 
interested in the m > this gives rise to a spacetime with positive Bondi 

mass. The spacetime is then determined by the evolution of P = P(u, on a 
background two-sphere. The evolution equation is often written as 

(hxP), u =--^-AK (2) 



This equation, the Robinson- Trautmau^ equation, is also known in the literature as 
the two-dimensional Calabi equation! 
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The steady state, corresponding to the Schwarzschild solution, is given by P = Pq = 
^ (a + b( + b( + cCC) , where ac — bb= 1. Note that the equilibrium value of P is not 
unique but includes a freedom corresponding to conformal motions on the sphere. 
This freedom was a contributing factor to the difficulty of proving existence of 
solutions of the Robinson- Trautman equation!] The condition ac—bb = 1 normalises 
the Gaussian curvature of S 2 to K = 1. 

The function P determines the induced 2-metric on each <S Ujr , which is given by 

2r 2 d(d( 2A 
9s = p2 = e 9o 

where go is the metric of S 2 , and e~ A = -j=p. It is possible to "factor out" the 
background S 2 geometry and write the Robinson- Trautman equation as 

e 2A A, u =-J-A ^ (3) 
12m 

where Ao is the Laplacian on S 2 . We thus solve the evolution equation on a back- 
ground sphere, say at r = 1. 
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2.2. Conserved quantities and Lyapunov functionals 

Several conserved integrals are implied by the evolution equations. First, we have 
the conservation of surface area of S u . r , and the conservation of the "irreducible 
mass" , as a consequence of this: 



A s = I u>i = e 2A uj = 47r 

/So„ JS 2 



If m I ,a 

Mi = - — / rnuii = — I e ujq = m. 
An ./.<,_. 47T ,/ S 2 



The Gaussian curvature, in terms of A, takes the form K = e~ 2A (l — AoA) which, 
making use of Stokes' theorem, immediately gives the conservation of the Euler 
number for S u y. 

XS = ^~ f Ke 2K Lo = 2. 
^ Js 2 

SingletonED gave an expression for the Bondi-Sachs mass of the spacetime: 

Mb = £L / e 3 % (4) 
4tt J S 2 

which is manifestly positive and in fact is bounded below by the irreducible mass, 
since by the Holder inequality 




That the Bondi mass is monotonically decreasing can be shown by differentiating 
Eq.(g) and using the Robinson- Trautman equation: 

i MB = ~h L eA ( P ° 2 ^ A )'f)'C ( P ° 2 (e- A ), c ),C < (5) 

These characteristics enabled Singletoni to show that the Bondi Mass is a Lya- 
punov functional for the Robinson- Trautman evolution, complementing the ear- 
lier work of Lukacs et. alcil who showed that the integral J s2 K 2 e 2A LJo is a Lya- 
punov functional for the Robinson- Trautman evolution. These -Lyapunov function- 
als played a key role in Chrusciel's semi-global existence proof.cl 
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3. The Past Apparent Horizon 



3.1. General Remarks 

Apparent horizons in black hole spacetimes have been the subject of much study in 
recent times. The motivation for the study of the future or past apparent horizons 
is that they provide a locally characterisable indication of the presence of an event 
horizon (future) or particle horizon (past). Indeed, in cases where the event horizon 
or particle horizon cannot be properly defined, the apparent horizon may be the 
only useful definition of the "surface" of the black hole or white hole.Q This arguably 
is the case in the Robinson- Trautman spacetimes, where J~ cannot be properly 
defined due.-to the instability of the Robinson- Trautman equation in the negative 
u direction.B Apparent horizons have been widely used ip_iuimerically generated 
spacetimes as an indicator of the presence of a black holejij'ta in lieu of being able 
to detect the event horizon in a local wayf It is believed that the presence of an 
apparent horizon indicates the presence of a nearby event horizon (where it exists) , 
and it has been shown thai—for all stationary black hole spacetimes the apparent 
and event horizons coincide. LJ Recent work has shown that the apparent horizon of a 
black hole has "thermodynamic" properties: in particular it has been shown that a 
future apparent horizon can only increase in area, while a past apparent horizon can 
only decrease in area, obeying a law like that of the first law of thermodynamics.^!!! 

The Robinson- Trautman spacetimes, as they settle down to the Schwarzschild solu- 
tion as u — > oo, would be expected to have a well behaved future apparent horizon 
coinciding with the event horizon. Due to the presence of outgoing radiation, we 
would expect the past apparent horizon to be more complex. As a general (and non 
rigorous) consideration, it can be seen quite clearly from the line element Eq.(|lf) 
that the curvature singularity at r — is a spacelike singularity i.e. that if the 
spacetime were complete in the past null direction (if J~ were complete) then 
there would exist a particle horizon - the past counterpart of the event horizon 
(see Fig.(P). Effectively, the past singularity being spacelike indicates that no null 
geodesies from J~ can reach the singularity, if the spacetime is to be causally well- 
behaved. Because the particle horizons is a globally defined object, it does not 
exist in Robinson- Trautman spacetimes. However, the spacelike nature of the past 
singularity suggests that there should always be a past apparent horizon in these 
spacetimes, which we believe to be causally well-behaved. 

3.2. Apparent Horizon Equations 

The concept of an apparent horizon rests on the definition of marginally trapped 
surface. We define an outer past marginally trapped surface, T~ , to be a spacelike 
two surface on which the ingoing future directed congruence of orthogonal null 
geodesies has vanishing divergence, while the outgoing future directed congruence 

* In recent work, Anninos et. aJ3 have been able to locate the event horizon in numerically 
generated spacetimes by integrating backwards from the stationary final state of the black hole 
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Figure 1: Penrose diagram for the Robinson- Trautman spacetime. 



of orthogonal null geodesies is diverging. We define the outer past apparent horizon, 
TiT , to be a hypersurface r = 5ft(it, C, C) such that its intersection with each u = uq 
slice is an outer marginally past trapped two-surface T u ~ , i.e. TL~ is foliated by the 
surfaces T u ~ . The past apparent horizon is illustrated in Fig.(Q). 

The equations describing the apparent horizon are constructed as follow. Consider 
a null tetrad (l a , n a , m a , fh a ) for the spacetime with l a tangent to a congruence of 
outgoing shearfree diverging null geodesies, and m a and fh a lying on a family of 
two surfaces %.: 



0, 
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P 2s ft, 



(6) 



where H is as in Eq.(Q) and r = 3?(u, £, Q defines the surfaces T u . The nonvanishing 
Newman-Penrose quantities are then given by: 



--, r = 7r = a + 
r 
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X" - -(m + P 2 5R, CC - ) + ^(2P 2 5ft, c 5ft, f ) 
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(7) 



From Eq.(§ and Eq.(0) it can be seen that the null vector l a is future directed, 
outward pointing, geodetic (k = 0), diverging (p = — -), null and orthogonal to the 
two-surface T u . That is, l a is tangent to the outgoing congruence of null geodesies 
normal to T u . The tetrad vector n a will be tangent at T u to the ingoing congruence 
of null geodesies normal to T u . The divergence of n a is given by Pe(/i) (= /i in 
this case). Thus the divergence of the ingoing congruence of null geodesies is given 
on T u by Jl, where the tilde is used here and subsequently to denote the restriction 
of an r-dependent quantity to the hypersurface r = 5ft(u, If A* = 0, 7^ are 

marginally past-trapped two-surfaces T u ~ and the hypersurface r = 5ft(u, £,C) is 
the past apparent horizon 7i~ . This condition gives us the first horizon equation: 



K 



2m 



AQnft) =0 



(8) 



This equation is attributed to PenroseEl. TodH examined it in more detail, proving 
that Eq.(||) has a unique C°° solution given a C°° background S u>r , and also that 
the surface T u ~ defined by the solution is in fact the outer boundary of past-trapped 
surfaces on u = Uq. 

The other equations are found by examining the embedding of the past apparent 
horizon in the spacetime; in particular by considering the normal vector to the 
apparent horizon. 

Let iV a be a one-form on TL~ defined by N a = — 5ft, u du + dr — 5ft,£ d( — 5ft, ^ dQ. Hence 
N a := g ab Ni, is a vector orthogonal to Tt~ and is given by 



N a 



n a - 5ft, u +H 



P 2 5ft, c 5R, c -~ 
jft 2 



I" 



P 2 5ft ? P 2s ft f 
d u - +2H)d r + —^d c + 
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5R 2 c 



(9) 
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Since the complex null vectors fh a and m are tangent to the two-sufaces which 
foliate the hypersurface 7i~ , the null vectors l a and n a are orthogonal to T u ~ . It 
can be seen that the vector 



Z" 



n a + \ $t,u +H 
d u + [ + 



P 2 5R, f 
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2P 2 K, c 5i,A , P 2 K, C - P 2 5R, C 



5ft 2 



SI? 2 



5ft 2 c 



(10) 



is orthogonal to N a and therefore is tangent to the hypersurface H . The "magni- 
tude" of iV a and Z a , from (|) and @, is given by 



-2 K,„ +P 



P 2 5ft, c 5ft, c -' 
5ft 2 



(11) 



Since from Eq. (||) /i = on each T u , the directional derivative of [i along the vector 
Z a tangent to Hr must vanish. (QQ) and (O) then imply 



Z a V a Jl = n a V a Jl - ^(N b N b )l a V a Jl = 



Substituting Eq. (^) into the Newman-Penrose equationsEl gives 



(12) 



n a V a /i = m a V a v + v(-t + a + 3/3) + nV - AA - 20 2 </> 2 



-0< 
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(13) 



(14) 



In (ph we have used v = —^(N a N a ) < and A = (^ i ), C . Combine @), @ 
and (|14|) to give 



1 f N N aS 

_Af^-)+2^(.V„.V") = -lAA 



To summarise, then, the apparent horizon in asymptotically flat Robinson- Trautman 
spacetimes can be described by the following equations: 



1m 

K- — =A(lnR) 



(15) 
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-K 



2m 



2P 2 5ft, f 



IP 



(16) 
(17) 



The i?e(/2) = equation Eq.(|l5|) is the primary equation required to locate the 
horizon. Tod's proof of existence and uniqueness of solutions for this equation 
indicates that there must be a marginally trapped two-surface on each slice u = uq. 
However, in order for TL~ , the union of these surfaces, to be considered an apparent 
horizon, we require that H~ be a non-timelike hypersurface, i.e. that N a N a is 
always non-negative. As shown by Collinsa and Haywardjj this is equivalent to 
the assumption that the directional derivative of the divergence of the ingoing null 
geodesies in the direction of future pointing outgoing null geodesies is negative. In 
fact, Hayward makes this assumption a part of the definition of an outer apparent 
horizon, and thus does not require a spacetime to be asymptotically flat in order for 
the outer apparent horizon to be selected. In the Robinson- Trautman spacetime, 
which is asymptotically flat, we can prove that this assumption is not necessary; 
using the maximum principle and Eq.(|l6|) it is straightforward to show that N a N a > 
(details of the proof will be given elsewhere). 

3.3. Properties of the horizon 

Equation Eq.([l5|) can be written in the form 

2me 3A = e <I, (l-A $) (18) 

where e* = Ke A . This gives a useful expression for the Gaussian curvature of the 
apparent horizon, Kt = e _2 *(l — Ao$). It is then straightforward to show that 
the marginally trapped surfaces have the topology of S 2 : 

Xt = ^- f K T e 2 *u = -!- / [1 - A $H = 2. (19) 

i-K Jg2 27T J S 2 

Furthermore, because = 2me 3A , Kt must be positive everywhere. Also, this 

relationship gives rise to some additional interesting integral quantities: 



me 2A LO = 47TTO <^> / if T e 3 * A ^o = 87rm 
S 2 Js 2 

me 3A ^wo = 2tt & [ K T e 2 *uj = 4tt 
s 2 Js 2 

me 3A u Q > / K T e 3<s> aj > 

s 2 Js 2 

4- I me\ < 4- f K T e 3<s, LO < (20) 

du J S 2 du J S 2 
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This is suggestive of the fact that the Robinson- Trautman equation can be rewritten 
entirely in terms of "horizon quantities". From Eq.([18|), we have 



■"-(s;' , < 1 -* # > 

and 

A A= ±A (<P + ln(l- A $)). 
Using the Robinson- Trautman equation in the form 

(e 2A ), W = --^-A (e- 2A [l-A„A]), 
12m 

we can construct an equation in $ alone. While rather messy, this formulation of the 
Robinson- Trautman equation is physically interesting as it describes in some sense 
the "evolution" of the apparent horizon of the white hole. The term "evolution" is 
used loosely here, since we expect the past apparent horizon to be non-timclikc. 

The properties of the area of the apparent horizon, At — J s2 e 2 *a->o, can a ls° be de- 
rived from the horizon equations, Eq.(|l5|) to Eq.jl^). Tod proved the isoperimetric 
inequality, 16nMg > At, for the area of the horizon.Ea In fact, 16ttM§ > At > 
167rm 2 , since, from the Holder Inequality and Eq. (po[) , 

1 2. 

A T * I — J > 4?r 
=>■ A T > 167rm 2 . (21) 



Since Eq.(|l6|) implies that the apparent horizon must be a non-timelike hypersur- 
face, i.e. that N a N a > always, we can show quite easily that the At must always 
decrease. From Eq.@ and Eq.(0), we have 

( e 2 *), u +A (e*- A ) = -e*+ A (N a N a ) (22) 

Thus 



d : A r = 4: I <' 1 '^ -/ e*" A |.V".V.,u, . 0. (23) 



du du 



Thus we have shown that the past apparent horizon TL~ is a non-timelike hyper- 
surface foliated by past marginally trapped surfaces 7~ u ~ whose surface area is a 
monotonically decreasing function of u. 
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4. Numerical modelling 



We are currently undertaking to solve the apparent horizon equations in conjunction 
with the Robinson- Trautman equation numerically, and to use the numerical model 
to investigate further some of the properties of the apparent horizon. Our results at 
this stage demonstrate the nice properties of the horizon equations. We are initially 
solving the axisymmetric equations: thus the system is 1+1 dimensional only. A 
full discussion of the numerical results will be presented elsewhere, but we give here 
a brief introduction to the work. 

In order to find the apparent horizon, we first solve the Robinson- Trautman equation 
to find the "background spacetime". For this, we follow the work of SinglctorO 
and Prager and Lunjij using exactly the Crank-Nicolson algorithm described of the 
latter. Next, on each 'time' slice u = uq, we solve Eq.(|l5|) using Newton-Raphson 
iteration. That is, for the values of P (or A) determined by the Robinson- Trautman 
equation, Eq.(|l5|) is solved to find the function 3? which describes the position of 
the marginally trapped surface. Because the convergence of the Newton-Raphson 
method is entirely dependent on making a good initial guess, we actually begin with 
the last timestep of the evolution, where the marginally trapped surface is known 
to be close to r — 2m, and solve for each time slice in turn, going backwards to 
u = 0, taking as an initial guess for each slice the solution on the previous slice. 

The example shown has the initial condition / = = l.+0.1Yi,o — 0.2Y2.0 - 0-3^3,0j 
where Y^o are spherical harmonics restricted to the axisymmetric case (i.e. Legendre 
polynomials). For convenience we set 12m = 1. Fig.(|])(a) shows a plot of / against 
x = cos 9 evolving through u. It can be seen that the system settles down to 
equilibrium {P = Pq) very quickly. Note that the equilibrium solution includes 
a component of the first harmonic in this case, so / is of the form / = a + bx. 
The function 3J, representing the position of the marginally trapped surface T u ~ 
is plotted in Fig.(Q)(b), also showing very smooth behaviour. 5R is not subject to 
this ambiguity in the final state, and settles down as it should to r — 2m, which 
corresponds to the Schwarzschild equilibrium. 

In Fig.(|J) the "horizon mass" , defined as Mj- — y jj?^, is plotted as a function of u 
with the Bondi mass and the irreducible mass, demonstrating the monotonicity of 
the horizon area, and the inequalities described above (Eq.(pi|)). 

5. Concluding Remarks 

In this paper we have demonstrated the existence of the past apparent horizon in the 
vacuum asymptotically flat Robinson- Trautman spacetimes, and described some pf 
its properties. Tod's proof of the existence of marginally past trapped surfaces,E3 
together with our result that they foliate a non-timelike hypersurface, allow us to 
call TL~ a past apparent horizon. It is pleasing to find that the Robinson- Trautman 
spacetimes, which are the simplest asymptotically flat radiating spacetimes, exhibit 
a well behaved apparent horizon structure. In more general radiating spacetimes, 
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Figure 3: Monotonic decay of mass quantities in the numerical evolution of the Robinson- 
Trautman spacetime. 

we would not expect the existence of apparent horizons to be guaranteed. This 
is demonstrated by examples such as the electrovac Robinson- TrautmanEiJ or the 
vacuum Bondi-Sachs spacetimes. 

We have also shown that the surface area of the past apparent horizon decreases 
monotonically with the retarded time u. This result is^ particular example of 
more general theorems about apparent horizon dynamicsm. However, we have not 
required all the assumptions used in the proof of those theorems. 

We have illustrated a method for describing apparent horizons analytically based on 
the construction of a suitable null tetrad for the spacetime - this method could be 
easily applied to other cases. In a future paper we will discuss further the numerical 
solution of the horizon equations, as a means of investigating the physics of these 
spacetimes. 
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